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Abstract 

Wireless networks with energy harvesting battery equipped nodes are quickly emerging as a viable 
option for future wireless networks with extended lifetime. Equally important to their counterpart in the 
design of energy harvesting radios are the design principles that this new networking paradigm calls 
for. In particular, unlike wireless networks considered to date, the energy replenishment process and 
the storage constraints of the rechargeable batteries need to be taken into account in designing efficient 
transmission strategies. In this work, such transmission policies for rechargeable nodes are considered, 
and optimum solutions for two related problems are identified. Specifically, the transmission policy that 
maximizes the short term throughput, i.e., the amount of data transmitted in a finite time horizon is 
found. In addition, the relation of this optimization problem to another, namely, the minimization of the 
transmission completion time for a given amount of data is demonstrated, which leads to the solution 
of the latter as well. The optimum transmission policies are identified under the constraints on energy 
causality, i.e., energy replenishment process, as well as the energy storage, i.e., battery capacity. For 
battery replenishment, a model with discrete packets of energy arrivals is considered. The necessary 
conditions that the throughput-optimal allocation satisfies are derived, and then the algorithm that finds 
the optimal transmission policy with respect to the short-term throughput and the minimum transmission 
completion time is given. Numerical results are presented to confirm the analytical findings. 



This work was supported by NSF Grant CNS 0964364. The conference paper [1], scheduled to appear at ICC 2011, contains 
technical material in part from the short term throughput optimization problem considered in this submission. 



I. Introduction 

With the wide deployment of battery powered wireless devices, prolonging the lifetime of 
wireless networks is becoming ever more critical [|2|. Systems powered with batteries suffer 
from a limited lifetime, whereas networks consisting of energy harvesting or rechargeable nodes 
can survive perpetually [3|. The performance that such systems can deliver is tied closely to 
efficient utilization of energy that is currently available, as well as that is to be harvested. Towards 
understanding fundamental performance limits with energy harvesting, in this work, we iden- 
tify optimal transmission policies for maximizing throughput and minimizing the transmission 
completion time of a wireless node with a replenishing energy source and a finite battery. 

There has been a significant amount of previous work in energy efficient communications and 
networking of nodes with non-rechargeable batteries, from the perspective of various protocol 
layers, or cross-layer approaches, e.g. [|4j-[[T2j. References relevant to this paper include Q, [[5J. 



Energy efficient deadline-constrained communications is considered in (4J resulting in a packet 
scheduling algorithm that achieves minimum energy. In pj, an insightful approach to energy- 
efficient rate control is developed and the problem with individual deadline or buffer constrained 
systems with data arrivals is solved. 

There has also been some recent interest in wireless networks with energy harvesting nodes. 



A queue stabilizing transmission policy is developed in [13 1 for a recharging battery powered 
transmitter. This is a modified adaptive backpressure policy that is shown to be asymptotically 
optimal for sufficiently large battery capacity. Reference [[141 calculates energy management 
policies that are throughput optimal or mean delay optimal, and shows that a greedy policy 
is optimal for both problems in low SNR regime. A discrete-time battery model is introduced 



in [15 1 for the problem of energy aware routing in wireless networks powered with renewable 
energy. Reference [ [T6| assumes a Markov model for battery state, and proposes a threshold 
algorithm to best utilize available energy to packets with different rewards. In |T7|, policies 
based on the energy-error probability tradeoff are developed to maximize successful transmission 
probability while minimizing probability of running out of energy in an energy harvesting body 



sensor network. For a solar powered wireless network, [18| lays out sleep/wake-up strategies 



for various factors and determines optimal parameters of the solar energy harvest based strategy 



using a bargaining game model. The most relevant reference work to the present paper is [19|, 
p0| , where the authors have considered an energy harvesting node, and found transmission 
policies that minimize the transmission completion time of a given amount of data. This work 
assumes an infinite capacity battery for the energy harvesting transmitter. 

In this paper, we consider the problem of maximizing the short-term throughput of an energy 
harvesting node, in other words, maximizing the data transferred under a deadline constraint. 



Unlike previous work [19|, we employ the realistic constraint that an energy harvesting battery 
must have finite energy storage capacity. We seek to find the optimum transmission policy under 
this energy storage constraint as well as the energy causality constraint. That is to say that, due 
to the finite battery at the transmitter, any received energy that overflows its capacity is lost, and, 



that energy can not be expended prior to being harvested. Like previous work [19|, we consider 
known energy arrivals, in order to find a bench mark solution for any transmission policy that 
considers communication under a deadline. We also show that the problem solved in this paper 
is closely related to the problem of transmission completion time minimization considered in 
p9| . Specifically, we show that the solution to the former is identical to that of the latter for 
the same parameters, providing a solution to the latter under battery storage constraints as well. 
The algorithm that finds the optimum power allocation to solve the former is developed. Then, 
using the relation between the two problems, a similar algorithm that finds the optimum power 
allocation to minimize the transmission completion time is presented. 

The remainder of the paper is organized as follows. The system model and the problem 
definition is presented in Section [TTJ. Section III solves the short-term throughput maximization 



problem and presents the optimal policy. Section IV addresses the completion time optimization 
problem. Numerical results are presented in Section [V] followed by concluding remarks in 
Section |VD 



II. System Model and Problem Definition 

We assume a single link continuous time system where a node transmits continuously and its 
rate can be varied at will via power control. Specifically, the transmitter can choose to transmit 
with power P{t) at any instant t, achieving a corresponding rate r(p(t)) where r(.) is a non- 
negative, increasing, strictly concave function that we will refer to as the power-rate function. 
Using a power-rate function of this form is fairly common (5J, [20|, and is clearly valid for 



the additive white Gaussian noise (AWGN) channel. Factors such as processing power, battery 
leakage or base energy consumption are not explicitly taken into account in this model, however, 
can be integrated in the power-rate function without violating its above stated properties. For 
example, a base energy consumption of power could be modeled as a shift in the power rate 
function horizontally by P b , or a processing power linear in transmission rate could be added to 
r(.) to yield a new power-rate function satisfying the same properties. 

We consider an energy harvesting system with a finite-capacity battery. That is to say that 
the battery can store up to E max . This storage capacity is considered to be constant throughout 
the transmission horizon, i.e., battery wear or fatigue is assumed to be much slower than the 
time scale of the problem. The property of keeping the battery energy below its capacity and 
non-negative, i.e. within [0, E max \, will be referred to as energy-feasibility in the sequel. 

We will consider the offline optimization problem. Hence the energy replenishment process 
is modeled as a discrete process with energy harvests of size E n arriving at time instances s n , 
as shown in Figure [T] where E n and s n take positive real values and are available non-causally 
to the transmitter at the beginning of transmission. The first energy arrival £"0 is conventionally 
at time s = 0, representing the initial energy in the battery. Due to the limited battery, if the 
harvested energy E n is larger than the available space in the battery at time s n , the battery is 
charged to maximum capacity and the remainder of the energy packet is discarded. Note that an 
instantaneous energy consumption requires infinite instantaneous power, which is not allowed 
by system definition. Thus, all harvested energy must first be stored before consumption. Based 
on this observation, it is safe to assume that E n < E max , as no more than an E max amount of 
a harvest can be utilized by the transmitter. Therefore, a truncation at E max for E n values will 



4 



be assumed in the sequel. 

An example to justify this model could be a solar powered surveillance network. A continuous 
operation would be necessary for surveillance purposes, and the nodes would realistically be 
operating on finite batteries. Once the nodes are deployed, the solar energy harvesting process 
would be predictable enough to have good estimates of E n and s n , yet significantly varying 
throughout the day. 

In this model, the amount of energy available at any time instant is constrained, either because 
a sufficient amount may not yet be harvested or cannot be stored in the battery. As such, we 
have an energy-feasibility constraint on the transmission policy. Specifically, a energy-feasible 
power allocation p(t) is a bounded non-negative function for the transmission power that ensures 
the battery state stays within [0,E max \. We express the set of energy-feasible power allocation 
functions as 



which is a convex set of functions, i.e., any convex combination of two energy feasible allocations 
is also energy feasible. Note that allocation functions which lead to a battery overflow are 
physically possible, but are discarded from the energy-feasible set due to being strictly suboptimal 
as shown in subsequent sections of this paper. 

Our goal is to find the optimum power allocation for maximizing the total data transferred 
from an energy harvesting node under a deadline constraint. The only constraint induced by the 
model on the power allocation is the feasibility set The objective is to maximize the total 
number of bits departed in the time interval [0, T] over p(t). Under a rate function r(p(t)) and 
a time deadline T the problem can be expressed as: 



n-l t > 



W = {p(t) \0<Y,E k - / p(t)dt 



< E, 



S n -1 <t'<S n } 



(1) 



PI : 



max 



p(t) 





(2) 



Next, we shall solve this problem. 
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III. Short-Term Throughput Maximization 

This section considers the throughput maximization problem for an energy harvesting node 
with infinite backlog and finite time. First, the necessary properties of the optimal policy are 
established. Then, an algorithm is presented to generate the policy that satisfies all the necessary 
conditions and the optimality of this policy is shown. 



A. Optimality conditions 

The following lemmas provide the necessary conditions for a power allocation policy to 
be optimal. These conditions provide valuable insight into the development of the algorithm 



proposed in Section III-B and the proof of its optimality. 



Lemma 1: Given the total amount of energy consumed in any time interval [ti, t 2 ), a constant- 
power transmission is throughput optimal, i.e., among all transmission policies consuming energy 
E, the constant-power transmission p(t) = j^- departs the largest number of bits. 

Proof: The proof is by contradiction. Assume that p(t) is non-uniform in (£1,^2) , i-e., 
for ti,t 2 G (ti,t 2 ), 7~i 7^ r 2 , we have p(ri) 7^ £>(r 2 ). Without loss of generality, assume that 
p(tx) < p(r 2 ). Keeping pit) unchanged at every other time instance, we define an alternative 
policy p'(t) as 

p(t)+5 te [ti,Ti + e] 

At) = \p(t)-6 te[T 2 ,r 2 + e} (3) 
pit) otherwise. 

where 5 > is an infinitesimal power displacement among two intervals of arbitrarily small 
duration e > 0. Clearly, the energy consumed by p(t) and p'(t) are equal. However calculating 
the information sent in [ti,t 2 ] yields 
-t 2 



B= r(p(t)) dt 



r(p(t))dt+ r(p(t))dt+ r(p(t)) dt + [r(p(rx)) + r(p(r 2 ))].e 

Jri+e J T2+e 
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fry fT2 fl 

< / r(p(t)) dt + r(p(t)) dt + r(p(t)) dt + (r(p(n) + 5) + r(p(r 2 ) - 5)).e 

(4) 

ta 

r (?/(*)) = £' 

L 

where the inequality in ([4]) follows from the definition of strict concavity for r(p). This shows 
that we can strictly improve the throughput by equalizing the two power values, and by extension 
uniformizing the power allocation. ■ 
Corollary 1: The optimal power allocation policy dictates the transmission power, and thus 
the rate, remains constant between energy arrivals. Consequently the power level can change 
only when a new energy packet arrives. This follows from Lemma [T] by considering the interval 
[s n , s n+ i] and introducing the feasibility definition in ([I]) to the optimization. However, given a 
total consumed energy in [s n ,s n+1 ], feasibility in the sense defined in ([I]) does not depend on 
the structure of p(t) in this interval. Thus the constant power policy is feasible, and therefore 
optimal for any given energy consumed. 

Lemma 2: Any power allocation policy yielding a battery overflow is strictly suboptimal. 

Proof: It is clear that the battery cannot possibly overflow without an energy arrival. Therefore, 
we start by assuming that the battery overflows at the energy arrival instant Sj. Let the power 
allocation allowing this overflow be p(t). Ei is less than or equal to E max by system model 
and causes an overflow, thus battery state at s~ is strictly positive. This implies that pit) can 
be increased by an infinitesimal amount 5 in (s^ — e, s») without violating energy-feasibility, 
which strictly increases the throughput. Since the excess energy required, 5.e, is recovered by 
the overflow at Sj, the remaining transmission schedule remains unchanged. Therefore, a power 
allocation that yields a battery overflow can not be optimal. ■ 

Lemma 3: In the optimal power allocation, the transmission power does not change unless 
the battery is either full or completely depleted. 
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Proof: Assume that at arbitrary time t transmission power changes, so that p(t~~) ^ p(t + ). 
Consider the interval [t — r,t + r], where the policy depletes a total energy of r.(p(t~) +p(t + )). 
Unless the battery is full or depleted at t, the energy feasibility constraints will be inactive in 
this interval. Let p* (t) = p ^ be the constant power transmission policy in [t — r, t + r] 

that expends the same total energy, which is also feasible for r sufficiently small, due to energy 
feasibility constraints being locally inactive around t. Then, p(t) can be replaced with p*(t) 
without altering the rest of the problem, and due to Lemma [T] departs strictly more bits. Thus, 
p(t) must be stay constant unless an energy feasibility constraint is active, i.e. the battery is 
either full or depleted. ■ 

Note that in the proof above, p*(t) might still be feasible when battery is full or depleted, 
specifically if it shifts the policy away from the active constraint. This observation forms the 
grounds of the next lemma: 

Lemma 4: For optimal power allocation, the change in power level p n at an energy arrival 
instant s n has to be nonnegative (nonpositive) if the battery is depleted (full) at that time instant. 

Proof: The proof is by contradiction. Consider the notation in the proof of Lemma [3} Battery 
being depleted at time t implies that the first inequality in ([!]) is active. However, if pitr) > p{t + ) 
holds, i.e. if transmission power is decreasing, then p*(t) < p(t). Hence replacing pit) with p*(t) 
on [t — r,t + r] increases the RHS of the active inequality, which is feasible. This shows that if 
the battery is depleted at t, a decreasing power policy at t can not be optimal. Similarly, battery 
being full at time t implies that the second inequality in ([T|) is active. This time, an increasing 
power policy, p(t~) < p(t + ), implies p*(t) > p(t); which renders replacing pit) with p*(t) on 
[t — r, t + t) feasible. This proves the second statement of the lemma. ■ 

Corollary 2: Lemmas [3] and [4] together imply that for optimal power allocation, transmission 
power decreases only at energy arrival instants when the battery is full and increases only at 
energy arrival instants when the battery is depleted. 

Lemma 5: The optimal power allocation expends all harvested energy by the end of the 
transmission. 



Proof: Assume that in an energy-feasible scenario, there is energy left after the transmission 
concludes. Without violating feasibility, this energy can be used in a nonzero time interval 
right before the end of transmission. This, due to power-rate function being increasing, strictly 
improves the throughput. Therefore the former policy can not be optimal. ■ 

B. Throughput Maximizing Policy 

First, we present the definitions of the variables used in the policy. From Lemmas [T]j4] in 
IH-A[ we know that the power allocation policy will consist of intervals of constant 



section 



power, changing at some energy harvest instances s n only. We denote the subsequence of s n at 
which the transmission power changes as i n . Let s nmax = T be a dummy harvest indicating the 
end of transmission. The energy of this dummy harvest is arbitrary as it cannot be used by the 
transmitter, and is set to E nmax = E max for convenience. The power allocation then has to be 
of the form 

(5) 

t > s nmax = T 

We remark that once the specifics of the interval t = [0, i\] is determined, the remainder of the 
problem can be considered as a separate throughput maximization problem. That is, given the 
duration of this interval and the amount of information sent in this epoch, ii.r(pi), it remains 
to solve for the optimal power allocation for a modified problem with energy arrival times shifted 
by ii, a new initial battery state ^2 k . s Kil E k — i\.p\, and a new deadline constraint of T — i\. 
This means that once the first time slot of the optimal allocation is identified, the remaining 
power levels can be found recursively with the same algorithm, using updated parameters and 
shifted arrival times. Thus, we shall focus on determining the optimal power level in the initial 
epoch. The modified optimization problem described above will be referred to as the shifted 
optimization problem. 

We define two sets of powers {p [l],p [2], ...} and {p max [l] } p max [2], ...}, where p [n) and 
Pmax[n] are the constant power levels that would result in an empty battery at s~ or a full 
battery at s+ respectively if employed in [0, s n ). Note that these levels need not be feasible, or 



even positive, but only serve comparison purposes. We then define the set P = {P[1],P[2], ...} 
with elements as the closed intervals P[n] = [p ma x[n],po[n}] between corresponding elements of 
the two sets p and p ma x- This translates to a range of constant power levels for the n th arrival 
that would be energy-feasible at s n when the feasibility at previous arrivals are disregarded. 
Thus, we have: 



V Ei. V n Ei - E 
Po[n\ = — — , Pmax[n\ = — — s n < T (6) 

S71 Sri. 



P[n] = [Pmax[n\,p [n}} = {p \ p max [n] < p < p [n}} s n < T (7) 

P[n ma x] = {po[n max \} s nmax = T (8) 

Here, (j8]) follows due to E nmax = E max yielding p [n max } = p max [n m ax]- Since Lemma [5] readily 
suggests that the energy should be depleted at s Umax , the convenient choice of E nmax is justified. 
Based on this definition of the feasible power range, it can be deduced that for a constant power 
transmission starting from t = to extend to the n th energy harvest without violating energy- 
feasibility, its power level should be contained in the range P[k] for k — 1, ...,n, i.e., the step 
should be feasible through all harvests it extends over. This yields an upper bound n ub on the 
length of the first constant power transmission, that can be calculated as 

n 

n ub = msix{n\ p| P[k] ^ 0, n = 1, 2, .., n max } (9) 

k=l 

as for the later harvests, a feasible constant power level contained in all previous feasible power 
ranges P[n) do not exist. This gives us a range for feasible constant power transmission levels 
and its maximum duration. 

Remark 1: The following observation explains the intuition behind the optimal algorithm 
that will be proposed. Given an energy harvesting scenario and the corresponding {poM} and 
{PmaxM} sets, assume that a constant transmission of power pi and duration s ni is feasible. 
This transmission then satisfies p\ G Dfeli an( ^ cannot extend beyond s riub , as it is rendered 
infeasible at s nub+ i by one of the constraints. However, which constraint causes this is an 
important indication of how the policy tends to change after s ni . 
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A constant transmission with power p 1 either over-depletes or overflows the battery at s nub+ i. 
The former case implies that the power level after s ni needs to decrease, and the latter implies 
that the power level needs to increase in order to conform to feasibility conditions at s nub . This 
can be verified by calculating updated values of Po[n] and p ma x[ n ) f° r a shifted problem after the 
first step of the policy is determined. Figure [2] shows how the boundaries of the interval P[n] 
for the shifted problem, shown in red, move away from the chosen transmission power pi when 
re-calculated starting from time s ni . This indicates that if an interval P[n] falls below or above 
a chosen transmission power, than the corresponding intervals of the shifted problem will again 
fall above and below respectively, forcing the power of the next step to increase and decrease 
respectively. By Corollary [2} an increase or decrease in power can occur only at a harvest point 
and with the battery depleted or full. Hence, the choice of p\ in the optimal policy is restricted 
to pofai] an d Pmax[ n i] respectively for the two cases in consideration. Observing the tendency 
of the power level to increase or decrease after s ni therefore determines the power level p\. 

The analysis and intuition presented so far aims at insights on the optimal power allocation. 
Based on these, we prove that the following algorithm determines the throughput maximizing 
power allocation for a transmitter node with initial energy E , energy arrivals E n at times s n , 
battery capacity E max and a deadline constraint T: 

Throughput Maximizing Algorithm, Al 

1) Find the upper bound on the length of the first constant transmission using|9j If n ub = n max , 
transmit with constant power (X}fc=cT Ek)/T until the end of transmission. 

2) Determine whether the next power interval P[n ub + 1] falls below or above n^o-^t^]- 

3) If P[n ub + 1] > ff k l b P[k], transmit with 

k = s ni , pi = PoM, where n x = max{n \p [n) E Q P[k]}, 

k=0 

If P[n ub + 1] < rifc=o P t^]' transmit with 

n 

ii = s ni , pi = PmaxW], where m = max{n \p max [n] E pj P[fc]}. 

fe=0 
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4) Repeat algorithm for the shifted problem with modified parameters 

E = Y1 Ek ~ il - Pl i T' = T - i h n' max = Umax - Tli, 
k=0 

s n — s n + ni — i±, E n = £7 n _|_ ni , n = 1, ... } n max . (10) 

A graphical description of the algorithm is provided in Figure [3] which shows the feasible 
energy tunnel of the energy harvesting transmitter. The upper solid boundary represents the 
cumulative energy harvested and the lower solid boundary is the upper boundary shifted down 
by an amount of E max . The cumulative energy spent by the power allocation algorithm forms 
a continuous line, and must stay within this tunnel to conform to energy feasibility. A power 
allocation that goes above the tunnel spends more energy than available, and one that goes below 
causes a battery overflow. Therefore, the set of energy-feasible power allocations ^ lies within 
this tunnel. Note that this graphical depiction is similar to that in [[5J, where a tunnel is formed 
for the data transferred. In contrast, here, we have an energy tunnel. 

In Figure [5J the sets p [n] and p max [n] correspond to the slopes of lines from the origin to 
each of the corner points in the tunnel as shown with dashed lines. The interval P[n] represents 
the slope set of lines passing through the n th opening in the tunnel, and is marked with an arc 
on the figure for n = 2. The first step of the algorithm determines the longest constant power 
transmission that stays within this tunnel. The second step determines whether the most distant 
point on a wall that a line through the origin can reach is an upper bound or a lower bound. 
This is accomplished by comparing the first unreachable opening with the last reachable one. 
Finally, the third step selects the longest feasible constant transmission that ends in one of the 
sets po[n) and p m(M [n], allowing a change in transmission power for the rest of the problem. 

Theorem 1: Algorithm Al yields the optimal power allocation policy. 

Proof: Proof is by contradiction. We begin by assuming that there exists an optimal power 
allocation step {pi},{«i} such that pi ^ pi or i\ < %\. There are four possible distinct cases 
over which we will show contradictions on optimality, labeled in Figure [4] as I- IV. 

Figure |4|a) shows a case where the algorithm chooses to transmit with power p\ = p max [ni] 
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for a duration i\ = s ni . This selection implies that point A lies on the upper tunnel boundary due 
to the outcome of the second step of Al. A power allocation policy can differ from this selection 
by either trespassing to the upper region I or the lower region II. A feasible power allocation 
forms a continuous line that stays within the tunnel throughout the transmission time. Therefore, 
a power allocation extending to region II must cross the algorithm's line again before s ni . At 
this crossing point, the algorithm's transmission departs strictly more bits than the alternative 
due to Lemma [Tj rendering any allocation extending to region II suboptimal. On the other hand, 
any policy extending to region I must cross the dashed line to be feasible, at which point a 
constant power transmission up to that point is feasible and performs strictly better. Therefore, 
any allocation extending to region I cannot be optimal. 

A parallel statement holds for the second case shown in Figure |4jb) where the algorithm 
chooses to transmit with power p 1 = Pot^i] instead. Similarly, due to the second step of algorithm, 
point B is known to lie on a lower boundary. A policy extending to region III has to cross the 
algorithms line at some point in order to complete a feasible path to the deadline. At this crossing 
point it proves itself suboptimal as the constant power transmission is strictly better. On the other 
hand a policy extending to region IV has to cross the dashed line, at which point a constant 
power transmission up to that point is both feasible and strictly better. Therefore no other policy 
can perform better than the constant power transmission shown with the solid line. 

The termination step follows an analogous logic. Due to Lemma [5} it is known that the 
policy must terminate with a completely depleted battery at deadline T. The termination step 
suggests a constant power transmission to this point should be employed whenever feasible. 
As the constant power transmission is the optimal solution for the unconstrained problem as in 
Lemma [TJ it emerges as the solution to the energy constrained problem as well, provided that 
it is feasible. Consequently any alternative power allocation departing from this constant power 
transmission would be suboptimal, as the two policies have to terminate at the same point at 
which the constant power transmission clearly departs more bits. 

In summary, any transmission policy that differs from the one found by the algorithm and 
still is energy feasible must be suboptimal. Consequently, our proposed algorithm's step has to 
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yield the optimal policy. ■ 
With the necessary conditions and algorithm Al derived for the throughput maximization 
problem, we move on to an alternative problem with similar results. 

IV. Transmission Completion Time Minimization Problem 

A previous problem studied with an energy harvesting transmitter is the transmission comple- 
tion time minimization problem, introduced in pOj. The transmission completion time minimiza- 



tion problem focuses on determining an energy feasible power allocation p(t) e that, given 
the total number of bits to send as B, finalizes the transmission in the shortest time possible. 
This problem has total transmission time T as an objective function to be minimized, and a 
throughput constraint as seen in the mathematical expression of the optimization problem: 

P2 : min T, s.t. B - [ r(p(t))dt < 0, p(t) E (11) 



o 



In contrast to its throughput maximization counterpart we considered in Section [HI] the time 
minimization problem appears to have a more complex form. Reference [20] solved this problem 



by employing lemmas similar to the lemmas in Section III- A , without the battery constraints. To 
address the model with finite battery capacity, an extension of the algorithm in p0| is possible. 
However, we will choose another direction that is simpler. 

Specifically, we note that although the throughput maximization problem differs significantly 
in structure from the transmission completion time minimization problem, their solutions are 
closely related. Theorem [2] states this relationship. 

Theorem 2: For a given energy harvesting scenario, the two optimization problems yield 
identical power allocation policies for matching time and bit constraints. In other words, if the 
maximum-throughput policy for time interval [0,T] departs a total of B bits, then the minimum- 
time policy for B bits completes the transmission at time T, and vice versa. 



Proof: First the Lagrangian dual problem of the time minimization problem in (11) is 



formulated in ( 12). The variables T and p{t) of the nested minimization problem are by definition 



independent. Keeping in mind that the solution to the dual problem satisfies u > 0, the inner 
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problem can be separated as in ( fT3] ). It can now be readily observed that the optimal power 
allocation function p*(t) arises from the solution of the maximization problem marked with PI 
for the optimal completion time T*. Observe that PI is identical to the throughput maximization 
problem ([2]). Therefore the solution of the completion time minimization problem is identical to 
the solution of the throughput maximization problem where the time constraint is the minimum 
transmission completion time T*. 

max ( min (T + u(B — [ r(p(t) )dt) j) (12) 

u>o Vp(t)e«p,VT Jo ' 

max ( min (T + u.B — u. max / r(p(t))dt )) (13) 

«>o V vt P(t)&pj ' ) 

v v ' 

PI 

■ 

Remark 2: An alternative proof is by making use of the monotonicity of the two problems 
in time and bits. The minimum-time problem yields a strictly larger completion time for more 
bits, while the maximum-throughput problem departs strictly more bits for a later deadline due 
to the strict concavity of the power-rate function. Assume that the maximum-throughput policy 
for interval [0, T] departs a total of B bits but the minimum-time policy for B bits at time 
completes the transmission at time T" ^ T. Consider the two cases T' > T and T' < T. The 
contradiction in the first case T" > T is trivial as the maximum-throughput policy departs the 
same number of bits in a shorter time; and thus the minimum-time policy cannot be optimal. In 
the second case, minimum-time policy achieves the same throughput in a shorter time T' < T 
indicating that there exists a shorter policy with the same throughput. Therefore strictly more than 
B bits can be sent in the larger time interval [0, T] and thus the suggested maximum-throughput 
policy cannot be optimal. 

Theorem [2] states that the solution to the former and latter problems are in fact identical. 
Therefore Lemmas [T]j5] that characterize the optimal power allocation also apply to the completion 
time minimization problem. We make use of this relationship to develop a modified algorithm 
that yields a throughput maximizing power allocation policy while departing exactly the desired 
number of bits, thus solving the completion time minimization problem. 
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Remark 3: The two algorithms, i.e. the one that yields the throughput maximizing allocation 
and the one that yields the completion time minimizing allocation, shall only differ in the 
termination condition. Throughout the time interval in which a constant power transmission 
until the end is not feasible, whether the end is defined by a deadline or end of a packet, the 
power allocation shall be identical. Consequently, the two algorithms work identically until the 
termination step is reached. 

The throughput maximization policy terminates at a certain time in contrast to the completion 
time minimization policy ending when a certain number of bits have been transmitted. Therefore 
the definition of the feasible transmission interval P[ra] in (|6])-([8]) needs to be modified as 

Y^ n ~ 1 rp rp _ zp 

Po [n] = ^^, Pmax [n] = ^ k= ° — A Po[n].s n <B (14) 

P[n] = {p\pmax[n] <p< p [n)}, Po[n].s n < B (15) 

P[n m ax] =Po[nmaz], Po[ n ]- s n ma x = B ( 16 ) 



where (16) suggests creating a virtual arrival point at s nmax that corresponds to the point for 
which a constant power transmission with the total energy harvested, ignoring energy constraints, 
transmits all bits. This is the adapted version of the virtual arrival point approach in ([8]), and 
can be interpreted as a candidate point for end of transmission that is selected if found feasible. 
The point s nmax lies between the two arrivals for which Po[n^\.s ni < B < Po[n i+ i\.s ni+1 holds, 
and can be found by solving the equation 

B = s n .r( ^ k=0 k ), s n . < s n (17) 

With the updated parameters we now provide the modified algorithm that yields the minimum 
transmission completion time, and prove that the algorithm yields the optimal solution. 

Transmission Completion Time Minimization Algorithm, A2 

1) Find the upper bound on the length of the first constant transmission using |9j If n ub = n max 
transmit with constant power (X]fc=cT Ek)/ Sn max until end of transmission. 

2) Determine whether the next power interval P[n u b + 1] falls below or above Plfc=o ^l^]- 
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3) If P[n ub + 1] > rifc=o P M' transmit with 

n 

ii = s ni , pi = PoM, where m = max{n \p Q [n] G Q P[k}}, 

k=0 

If P[n ub + 1] < Dfc=o p [ fc ]' transmit with 

n 

h = s ni , pi= Pmax[ni], where m = max{n \p max [n] E f~}P{k]}. 

k=0 

4) Repeat algorithm for the shifted problem with modified parameters 

ni 

K = Ek ~~ B' = B — r(pi).s ni} rt max = n max -i x , 

k=0 

S n = s n+ni ^1) -B^ = -^n+ni > Ti = 1, ...,Tl max . (18) 

Theorem 3: Algorithm A2 gives the optimal power allocation scheme for the transmission 
completion time minimization problem. 

Proof: We shall make use of the relation between the two problems to simplify this proof. If 
the suggested power allocation with completion time T* is identical to the throughput maximizing 
allocation with a time constraint T*, we can state that no other allocation with time constraint 
less than T* can depart B bits, and therefore the allocation is optimal in time minimization. 

This algorithm appears to differ from the throughput maximizing algorithm presented in 
Section [HI] at only two points: (|T4|)-([T6|) and (18). However, ( fTS] ) only affects the termination 



condition through n max in ( 16). Consequently the two allocations found using the two algorithms 
are identical until either one reaches a termination step. The termination step of the time 
minimization algorithm is reached when there exists a feasible constant power transmission step 
that departs all the remaining bits by T*. The presence of this step implies that the matching 
throughput maximization problem with deadline T* would at the same time have the same 
feasible constant power transmission opportunity and terminate. Conversely if there does not 
exist a feasible last step departing all remaining bits, then T* is unreachable by a constant power 
transmission and the throughput maximization problem also fails to terminate. Thus we can state 
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that the matching algorithms terminating at the same time T* would reach the termination step 
at the same time. 

The termination step on the other hand suggests a constant power transmission until the end 
of transmission at T*. As the previous power policies for the two algorithms are identical, the 
energies left for the last step are equal. Therefore for any termination step decided by the time 
minimizing algorithm ending at T*, the throughput maximizing policy with the deadline T* 
is forced to take the same step mainly because a constant power transmission is feasible with 
the available energy. As a result, this algorithm yields a power allocation that is identical to a 
maximum throughput solution, and is therefore optimal by Theorem [2j ■ 



V. Simulation Results 

In this section, we present simulation results to demonstrate the behavior and the performance 
of the algorithms. First, we have a sample simulation run of the algorithm presented in Sec- 
tion 



III-B The simulation is performed for a node with E max = 10 units with energy arrivals of 
E n = {2, 1, 6, 4, 8, 1} units at times s n = {0, 2, 4, 5, 7, 11} respectively. The energy arrivals are 
shown in Figure [5^ a). The energy tunnel of this system is shown in Figure [5jb) along with the 
result of the throughput maximization algorithm for a deadline constraint T = 12. 

At the first step of the algorithm, the sets p [n] and r max [n] are calculated and the corresponding 
intervals P[n] are determined. These intervals are displayed in Figure [5jc) with vertical lines. 
The upper bound on n\ is determined to be n u b = 3 as the third interval P[4] falls outside 
P[l] n P[2] fl P[3]. Based on the position of P[4], the longest feasible transmission with power 
in po[n] is picked, and the first step of the algorithm is decided as p\ with length s ni . When new 
intervals for the shifted problem are calculated beyond s ni , an opposite position for P[n ub + 1] 
is observed and a transmission within p m ax[n] is picked instead. Finally, for the next shifted 
problem, the constant power transmission is found to be feasible and the algorithm is terminated. 

Observe that no other power allocation policy within the energy feasible tunnel can perform 
better than this. This can be proved by stating that any different power curve reaching from the 
origin to the end point at T = 12 must cross this line at least once, at which point it has spent 



18 



the same amount of energy while departing strictly less bits. 

Assuming the same energy arrival scenario, we next observe the performance of the algorithm 



presented in Section III-B and the algorithm presented in Section IV for a range of deadline and 
packet size parameters in Figure [6] As predicted by Theorem [2| the curves match up perfectly, 
verifying the relationship between the two problems. A point (Bi,T{) on this curve corresponds 
to a power allocation that solves both of the problems in consideration simultaneously: The 
throughput maximization problem for a given deadline T t is solved by departing a maximum of 
B\ bits, while the transmission completion time minimization problem for a given packet size of 
Bi bits is solved by completing the transmission in T\ with an identical power allocation policy. 
Another noteworthy observation from Figure [6] is the strictly increasing nature of both of the 
problems with respect to their parameters. As a consequence, more time is required to depart a 
longer packet, or similarly more bits can be transmitted when given a more lenient deadline. 

After observing the behavior of the algorithms on a smaller time scale, we simulate for 
longer realizations in order to evaluate the average long term performance of the algorithms. 
Assuming a battery capacity E max = 100 we generate energy arrivals randomly with energy 
packet size distributed uniformly in [0, E max ] and inter-arrival times distributed exponentially 
with an average of 5 seconds. In this setting, Figure [7] compares our optimal power allocation 
algorithm for the throughput maximization problem with two alternative power allocations for 
a deadline of T = 10000 sec. First corresponds to transmission without regard to any battery 
or arrival constraints. This is the performance of a "traditional" transmitter with no energy 
harvesting, and is presented as an upper bound for our model where we are bound to conform to 
energy feasibility and battery constraints. Second is an alternative algorithmic approach named 
the on-off algorithm. It is based on the fact that for a strictly concave power-rate relationship, 
constant power transmission is the most efficient. The policy is such that, the transmitter operates 
with constant power when energy is available, and shuts off when energy is depleted. The 
constant power level is determined from the average energy arrival rate. As seen in Figure [7J the 
energy constraints of the problem result in a performance loss with respect to having perpetual 
energy. However, a major portion of this loss can be recovered with the incorporation of arrival 
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information and use of variable power transmission. The optimal offline algorithm we present 
provides this remedy, and significantly outperforms the greedy on-off algorithm. 

Throughout further simulation runs involving various battery capacity and arrival statistics 
within a larger number of randomly generated arrivals, we observed that the optimal offline 
algorithms provided improvements over simpler policies such as the on-off transmission algo- 
rithm suggested above. The improvements were especially notable for the cases of relatively 
small battery capacity E max , large energy arrival amounts E n and large standard deviation on 
the arrival process. Therefore, it can be stated that employing optimal power allocation is more 
beneficial for systems with limited battery capacity or energy arrivals with large variations. 

VI. Conclusion 

In this paper, we have solved the short-term throughput maximization problem for a link 
with an energy harvesting transmitter, and limited energy storage capacity. Furthermore, we 
have shown that obtaining the maximum amount of data transferred by a given deadline, is also 
equivalent to solving for the minimum completion time (the deadline) given this amount of data. 
We have proposed the algorithms that yield the optimal solution of both problems and proved 
their optimality. 

The findings of this paper provide insight to developing optimal transmission policies for nodes 
that have some notion of when and how much energy they can harvest. Future directions include 
developing online power allocation algorithms for causal systems or systems with stochastic 
future energy harvests. There are also plenty of possible extensions for multiterminal system 
models from the single link model considered in this paper. 
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Fig. 1: Energy harvesting model. 
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Fig. 2: Pictorial depiction of the behavior of feasible power ranges for future arrivals. 
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Fig. 4: The optimal first step and suboptimal regions. 
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Fig. 5: Arrival scenario and simulation results for the sample run. 
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Fig. 6: Overlay ed plot of throughput by deadline and completion time by packet size for the 
energy harvest scenario in Figure |5j 
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Fig. 7: Performance comparison of transmission policies for a transmission of T = 10000s with 
E m ax = 100, E n distributed uniformly in [0,E max ], and interarrival times s n — s n -\ distributed 
exponentially with a mean of 5 seconds. 



